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' Abstract 

In this article is proved the existence and uniqueness of a minimizer of 
the energy for the non-relativistic one electron Pauli-Fierz model, within the 
class of pure quasifree states. The minimum of the energy on pure quasifree 
states coincides with the minimum of the energy on quasifree states. Infrared 
and ultraviolet cutoffs are assumed, along with sufficiently small coupling 
<^ I constant and momentum of the dressed electron. A perturbative expression 

\0 • of the minimum of the energy on quasifree states for a small momentum 

^ , of the dressed electron and small coupling constant is then given. We also 

^ I express the Lagrange equation for the minimizer, in terms of the generalized 

one particle density matrix of the pure quasifree state. 

o ■ 
m 

I Introduction 



X ■ I.l The Hamiltonian 



According to the Standard Model ofNonrelativistic Quantum Electrodynamics [|2l 
the unitary time evolution of a free nonrelativistic particle coupled to the quantized 
radiation field is generated by the Hamiltonian 

H, := l(lV.-A(f + (LI) 
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acting on the Hilbert space Iv^(M^; ^) of square-integrable functions with values 
in the photon Fock space 

oo 

d := ■■= 05-?^ (2), (1.2) 

n=0 

where (Z) = C ■ i7 is the vacuum sector and the n-photon sector 5+ ^ (2) = 
S{Z®'^) is the subspace of totally symmetric vectors on the ri-fold tensor product 
of the one-photon Hilbert space 

Z = [fe L'^iS^x^ C ® M.^) I yk G 5^,A a.e. : k ■ f{k) = O} (1.3) 

of square-integrable, transversal vector fields which are supported in the momen- 
tum shell 

Sa,A ■■= {keM.^\a< \k\ < A}, (1.4) 

where 0<o'<A<ooare infrared and ultraviolet cutoffs, respectively, reflecting 
our choice of gauge, namely, the Coulomb gauge. It is convenient to fix real 
polarization vectors €±{k) E M.^ such that {e+{k),e^(k), 7|:} C form a right- 
handed orthonormal basis (Dreibein) and replace (Oi l by 

Z = L\S,^^ X Z2), (1.5) 

with the understanding that f{k) = f{k, +) + e_ f{k, —). 
In (IlTI) the energy of the photon field is represented by 

Hf = J \k\a*{k)a{k)dk, (1.6) 

where / f{k)dk := Er=± /s, ^ /(^' ^) {«(^), «*(^)}fce5<,,AxZ2 are the 

usual boson creation and annihilation operators constituing a Fock representation 
of the CCR on^, i.e., 

[a{k),a{k')] =[a*{k),a*{k')] = 0, (1.7) 
[a{k) , a*{k')] =6{k-k')l, a{k)n = 0, (1.8) 

for all k, k' G ^^-^a x ^2- The magnetic vector potential A(x) is given by 
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with k = (k, t) eR^ X Z2, 

G{k,T) := ger{k)\k\-'/\ (I.IO) 

and g E M. being the coupling constant. In our units, the mass of the particle and 
the speed of light equal one, so the coupling constant is given as g = -^y/a, with 
a ~ 1/137 being Sommerfeld's fine structure constant. 

The Hamiltonian Hg preserves (i.e., commutes with) the total momentum op- 
erator p = 4 Vx • + -P/ of the system, where 

Pf = J ka*{k)a{k)dk (1. 11) 

is the photon field momentum. This fact allows us to eliminate the particle degree 
of freedom. More specifically, introducing the unitary 



[x 



(27r)3/2' 
(1.12) 



one finds that 



where 



iJHgW = / Hg^^d^P, (1.13) 



Hg,p = l{Pf + HO) -P)' + Hf (1. 14) 

is a selfadjoint operator on dom(ifQ the natural domain of i/g g = |ij + Hf. 

1.2 Ground State Energy 

Due to (|I.13I) . all spectral properties of Hg are obtained from those of {Hg^^}^^]g^3. 
Of particular physical interest is the mass shell for fixed total momentum p E M^, 
coupling constant g > 0, and infrared and ultraviolet cutoffs < a < A < 00, 
i.e., the value of the ground state energy 

E^,{g,p,cr,A) := mia[Hg^j,] > (1.15) 

and the corresponding ground states (or approximate ground states). 

We express the ground state energy in terms of density matrices with finite 
energy expectation value and accordingly introduce 

^ := [pE C\d) I P > 0, Trj[p] = 1, pH,^, H,^p E C\d)}, (1.16) 
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so that the Rayleigh-Ritz principle appears in the form 



(1.17) 



Note that TY^[pHg^p] = Tt^Ip^"'^ Hg^^^p^ for all < /3 < 1, due to our assump- 
tion pH^ g,H^^p eC\d)- 

It is not difficult to see that the ground state energy is already obtained as an 
infimum over all density matrices 



Dm := <pe 



pNf,Nfp eC\d) 



of finite photon number expectation value, where 



Nf 



a*{k) a{k) dk 



is the photon number operator. Indeed, if cr > then 



H„ri > Hf > aN 



(1.18) 



(1.19) 



(1.20) 



and 2)971 = 2)971 is automatic. Furthermore, if cr = then it is not hard to see [13 
that Egs{g, p, 0, A) = limo-\^o Eg^ig, P, cr, A), by using the standard relative bound 



'^<A6)ij\\ < O(0r)||(ify,<. + 1)1/2^ 



(1.21) 



where A<o^(0) and -f^/,<cr are the quantized magnetic vector potential and field 
energy, respectively, for momenta below a. So, for all < cr < A < oo, we have 
that 

Eg,(g,p,a,A) = inf {xr^ [p i/,,p^(cr, A)] | p G 2)97t}, (1.22) 

indeed. If the infimum (11.221) is attained at pgs{g,p, a, A) E DTI then we call 
Pgs{g, p, cr, A) a ground state of Hg^^{a, A). 

Since DDJl is convex, we may restrict the density matrices in (11.221 ) to vary 
only over pure density matrices, 

^gs {g,p,(r,A) = inf | Tr ^ [p Hg^p{a, A 

where pure density matrices are those of rank one. 



pD97t := {pe D9JI 



P e pDTl 

3^ G ll^ll = 1 : P= |^)(^ 



and 



(1.23) 
(1.24) 

pD97t := DTI n pD97t. (1.25) 

Another class of states that play an important role in our work is the set of centered 
density matrices, 

cDTl := ^p e DTI VfeZ: Tr^ [p a*(/)] = o}. (1.26) 
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1.3 Bogolubov-Hartree-Fock Energy 

The determination of Egs{g,p) and the corresponding ground state Pgsig,p) G 
'DVyt (provided the infimum is attained) is a difficult task. In this paper we rather 
study approximations to Egs{g,p) and Pgs{g-,p) that we borrow from the quantum 
mechanics of atoms and molecules, namely, the Bogolubov-Hartree-Fock (BHF) 
approximation. We define the BHF energy as 

EBHF{g,P,<y,^) = inf |Tr5[piJ,,p-(a,A)] | p E 0^}, (1.27) 

with corresponding BHF ground state(s) pBHpig,?, A) E £l^, determined by 

Tr^[pBHFig,P,(T,A)Hg./a,A)] = EBHFig,P,(T,A), (1.28) 

where 

0^ := jpeDOJl pisquasifreej C 2)971 (1.29) 

denotes the subset of quasifree density matrices. A density matrix p E DTI 
is called quasifree, if there exist fp E Z and a positive, self-adjoint operator 
hp = h*^ > on Z such that 

{WiV2f/t))^ := Tr^[pWiV2f/i)] = exp [2tlm{fp\f) - + hp)f) , 

(1.30) 

for all f E Z, where 

Wif) := exp[t^f)] := exp [^(a*(/) + a(/))] (1.31) 

denotes the Weyl operator corresponding to / and we write expectation values 
w.r.t. the density matrix p as {■)p. 

There are several important facts about quasifree density matrices, which do 
not hold true for general density matrices in DTI. See, e.g., [l3l|9l|5l[6l. The first 
such fact is that if p G is a quasifree density matrix then so is W{g)*pW{g) E 
for any g E Z, as follows from the Weyl commutation relations 

yf.gEZ: Wif) Wig) = e-^^'^^f\^^Wif + g). (1.32) 

Choosing g := -i\plfp, we find that Wi-iy/2fp)*pWi-iy/2fp) is a centered 
quasifree density matrix, i.e., 

WiV2fp/i)* pWiV2fp/i) E c£ld := mncDTl. (1.33) 

Next, we formulate a characterization of centered quasifree density matrices. 
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Lemma I.l. Let p G c2)9Jl be a centered density matrix and denote {A)p : = 
TT^{pA}. Then (i) {ii) ^ {in), where 

(i) p E til's is centered and quasifree; 

(ii) All odd correlation functions and all even trunctated correlation functions 
of p vanish, i.e., for all N G 'Nandipi, . . . , ip2N G Z, let either hn := a*{ipn) 
or hn := a{(pn),for all 1 < n < 2N. Then (61 ■ • • &2Af-i)p = and 

{hlh2- ■ ■h2N) p = 5Z (^'r(l) &7r(2))^ ■ • ■ (&7r(27V-l) &7r(27V))p, (1-34) 

where '^2N denotes the set of pairings, i.e., the set of all permutations vr G 
&2N of2N elements such that 7r(2n — 1) < 7r(2n + 1) and 7r(2n — 1) < 
n{2n),for all 1 < n < N — 1 and 1 < n < N, respectively. 

( Hi) There exist two commuting quadratic, semibounded Hamiltonians 

(1.35) 
(1.36) 



with B = B* > 0, C = C'^ G C?{Z), where {V'ijieN ^ Z is an or- 
thonormal basis, such that exp{—H — f3H') is trace class, for all (5 < oo, 
and 

_ f Tr,[A eM-H-PH')] \ 
^^^^ - rooX Tr,[eM-H-m] J' ^ ^ 

for all A e Bid). 

Eq. (11.331) and the vanishing (ii) of the truncated correlation functions of a 
centered quasifree state imply that any quasifree state p G 05 is completely de- 
termined by its one-point function {a{Lp))p and its two-point function (one-particle 
reduced density matrix) 

where the operators 7^, G B{Z) are defined as 

7p^) := (a*(^)a(<^))p and {if, a pip) := {a{y^) a{Jip)) p, (1.39) 
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and i7 : Z — 7- Z is a conjugation. The positivity of the density matrix p im- 
plies that r[7p, Qfp] > and, in particular, 7p > 0, too. Moreover, the additional 
finiteness of the particle number expectation value, which distinguishes DDJl from 
DDJl, ensures that 7^ G C^iZ) is trace-class, namely, 

Tr^bp] = {Nf)p < 00, (1.40) 

and that «p G C^i^Z) is Hilbert-Schmidt. 

Similar to (|I.24I) - (|I.25I) . we introduce pure quasifree density matrices, 

pm ■■= iid n pmn. (i.4i) 

A subset of of special interest is given by coherent states, i.e., pure quasifree 
states of the form \W{-i^f)n){W{-i^f)n\, which we collect in 

coh := {\W{-iV2f)n){W{-iV2f)Vt\ \f eZ]. (1.42) 

For these, 7p = = 0. 

Conversely, if 7 G C\{Z) is a positive trace-class operator and a G C^iZ) 
is a Hilbert-Schmidt operator such that r[7, a] > is positive then there exists a 
unique centered quasifree density matrix p G tO."^ such that 7 = 7p and a = ap 
are its one-particle reduced density matrices. 

Summarizing these two relations, the set O."^ of quasifree density matrices is 
in one-to-one correspondence to the convex set 

1-pdm := |(/,7,a) G Z©£^(Z) ©£2(2)1 r[7,a]>o}. (1.43) 

Note that coherent states correspond to elements of 1 — pdm of the form (/, 0, 0). 

Next, we observe in accordance with (11.431) that, if p G is quasifree then 
its energy expectation value {Hg^^)p is a functional of (/p, 7p, dtp), namely, 

{Hg,p)p = Sgdfp^P^^p): (1-44) 

where 



7, «) 
4{(T 

+ Tr[7A; ■ 7A;] + Tr[a*A; ■ ak] + TT[\k\'^-f] 



i{(Tr[7^] + rkf + 2Re(rG) - p)'' 



+ 2Re((G + kf) a{G + kf)) + Tr[(27 + 1)(G + kf) ■ (G + A;/)*]} 
+ Tr[7|fc|] + r|^|/. (1.45) 

Now we are in position to formulate our main results. 
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Theorem 1.2. Let 0<cr<A< oo, and p G M^, |p| < 1. Minimizing 

the energy over quasifree states is the same as minimizing the energy over pure 
quasifree states, i.e., 

EBHF{g,P,(r,-^) ■= ini TilHg^j^p] = inf Tr[iJg,p-p] . (1.46) 

Theorem 1.3 (Coherent States Case). There exists a universal constant C < 00 
such that, for < a < A < 00, ln(A + 2) < C and \p\ < 1/3, there exists 
a unique fg^^ which minimizes Sg^^{f) = £g^^{f, 0,0) in L'^{S„^\ x Z2, {\\k\'^ + 
\k\)dk). 

1. The minimizer fg^^ solves the system of equations 

fg,p 



^ %,p =P- 2Re(/* .G) - fl^kfg^j; , 
with \ug^^\ < \p\. 

2. ForO < a < A< 00, 

, £g,p{f)= inf £g,pif) = £a,pifg,p) , 

and for < cr < A < 00, fg^f G L'^{Sa,A x Z2). 

3. For fixed g, a. A, and small values of \p\, we have that 

£gAf9,p) = ^aA^) - P- G*-^^ ;4 ■ p + ) . 

^\k\^ + \k\ + 2G ■ G* 

We summarize in Theorem lL4l the information obtained in SectionslVllto lVIIIl 

Theorem 1.4 (Quasifree States Case). Let < a < A < 00. There exists C > 
(possibly depending on a and A) such that for all\g\,\p\ < C, there exists a 
unique {fg^p, 'Jg^p, (ig^p) which minimizes the energy SgAf^ I1 

1. The dependence of {fg,-^, ■jg^f, dig^ on {g, p) is smooth. 

2. The functions {fg,^, 7g,p', Cig,p) satisfy 

fg,p={l\k\' + \k\)-'p.G + 0{\\{g,p}r), 
ag,,^=-S^\G-G*) + 0{\\{g,m'), 

7g,p "I" 7g,p ~ ^g,p^g,p ' 

where S acts on the kernel KA{k, /c') of a Hilbert-Schmidt operator A as the 
multiplication by A; ■ k' + + \k\ + + \k'\. 
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3. For fixed a, A, and small values of \g\ and we have that 

= £g,,iO, 0, 0) - g'\p\'C2,2{'y, A) - g'C,,o{a, A) + 0{\\ {g, f ) , 
as {g,p) 0, with C2,2(o-, A) = {27r'^ - |7r) In(^) and C4,o(cr, A) > 0. 

4. The minimizer {fg,^, 'Jg^f, cng^p) satisfies (we drop the g, p indexes to simplify 
the notation) 

M(7, u)f = -ikl + \l)-u)-G-k-a{G + kf), 

A{X)a = -(G + kf) ■ {G + kfY, 
7 + 7^ = a a* , 

POO 

A:= / e-'^'^+^\M{-f,u) + {G + kf)-{G + kfy)e-'^-2+^Ut, 
Jo 

u:=p-TT[^k]-f*kf-2Re{f*G), 

with 

M{'y,u) := -\k\'^ + \k\ — k ■ u + k ■ ■jk, 
A{\)a := ka ■ k + Xce + ceX. 

Remark 1.5. In the coherent states case the formula 

£aAf9,p) = ^.,p(0) - /bf A) + 0{\\{g,p)\\') , 

holds and can easily be compared to the quasifree state case. 

Remark 1.6. Although Theorem |L4]is formulated in terms of the one-particle re- 
duced density matrix Tp and its constituents 7p and dtp, it turns out to be more con- 
venient to parametrize the pureness constraint 7p + 7^ = apCe*p in terms of an anti- 
linear Hilbert-Schmidt operator f which is chosen such that 7p = i(cosh(2f) — 1), 
ftp = |sinh(2f) J, where J : f E L^iS^^A x Z2) ^ / G L\S^^a x Z2). This is 
explained in detail in Section Hill 

Outline of the article We introduce our notation to describe the second quan- 
tization framework in Section HIl Section [III] introduces two parametrizations of 
pure quasifree states and contains the proof of Theorem |L2] The energy func- 
tional for a fixed value of the momentum p of the dressed electron is computed 
in Section |IVl along with some positivity properties of the different parts of the 
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energy. From Section |V] on we tacitly assume that the coupling constant > 
is small. The energy is then minimized in the particular case of coherent states 
in Section |Vl providing a first upper bound to the energy of the ground state and 
a proof of Theorem |L3] The existence and uniqueness of a minimizer among 
the class of pure quasifree state is then proven in Section IVT] provided \p\ is small 
enough. The first terms of a perturbative expansion for small g and p of the energy 
at the minimizer is computed in Section IVIII Finally the Lagrange equations as- 
sociated with the problem of minimization in the generalized one particle density 
matrix variables are presented in Section IVIIII 



II Second Quantization 

In this section Z denotes a C-Hilbert space with a scalar product C-linear in the 
right variable and C-antilinear in the left variable. 

Let B{X; Y) be the space of bounded operators between two Banach spaces X 
and Y, and C^{Z) the space of trace class operators on Z. Given two C-Hilbert 
spaces {Zj, {■, j = 1,2 and a bounded linear operator A : Zi Z2, set 
A* : Z2 ^ Z\ to be the operator such that 



V2;i G Zi, Z2 e Z2, {Z2, Azi)2 = {Zi, A*Z2)l , 

and ReA := |(A © A*), ImA := ^(A © {-A*)) G B{Zi, Z2) © B{Z2, Zi). 
Example 11.1. For z, z' G Z, 

{z,z) = z z . 
The adjoint of a bounded operator Aon Z h A* . 

The symmetrization operator 5„ on Z®" is the orthogonal projection defined 
by ^ 

Sn{Zi ®---®Zn) = —. y^Z^^®---®Z^^ 



ni 



and extension by linearity and continuity. The symmetric tensor product for vec- 
tors is 2i V ^2 = 'Sni+n2izi ® ^2) ^ud morc generally for operators is AiW A2 = 
Sg,+g, o {Ai © A2) o Sp,+p, for Aj G B{Z^Pi; Z^"^). We set 

:= gP'^ := 

Definition II.2. The symmetric Fock space on a Hilbert space Z is defined to be 

00 

Vn 



n=0 
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where := CVl, Vl being the normalized vacuum vector. 

For a linear operator C on Z such that ||C||b(2) < 1, let r(C) defined on 
each by C^" and extended by continuity to the symmetric Fock space on Z. 

For an operator X on Z, the second quantization dT(X) of X is defined on 
each Z^" by 

dT{X) 



= n 1^"-^ V X 



and extended by linearity to 0|^>o The number operator is Nf = dT{lz). 
For a vector / in the creation and annihilation operators in / are the linear 
operators such that a{f)^l = 0, a*(f)n = f, and 

= Mr 9) ^7^""' , and a*(/)^/^" = V ^7^" , (11.47) 

for all g E Z.By the polarization identity 

Eq. (III.47I) extends to and hence also to 0">o 2^"- They satisfy the canoni- 
cal commutation relations [a{f),a*{g)] = f*g, [a{f),a{g)] = [a* (f) , a* (g)] = 0. 

The self-adjoint field operator associated to / is $(/) = :^('^*(/) + '^(/))- 
For more details on the second quantization see the book of Berezin [4] . 

A dot "■" denotes an operation analogous to the scalar product in R^. For every 
two objects a = (ai, a2, as) and b = 62, with three components such that 
the products ajbj are well defined 

3 

a - h := ajhj . 



Example IL3. Withp e GeZ'\ke {B^'^' 



3 





3 


k ■ p = 


3 


P- 


G 




















3 










k-G = 




G* 


■ k 
















3 




3 






G* 


= J2g,g*^b''\ 


G*-G = 


J2G*G,eC, 












3=1 







3 
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where for an object with three components a — (oi, 02, 03) such that a* is well- 
defined, a* := (a*, 03, Og). We sometimes use the notation p'^ = or k '"^ — 

\k\\ 

And with another product, such as the symmetric tensor product V, 

Recall that the Weyl operators are the unitary operators W{f) — exp(i$(/)) 
satisfying the relations 

Vzi, Z2 e 2^ : W{z{)W{z2) = e-i^'^^'*'''^W{zi + z^) , (11.48) 

2 ^ Vn 

e Z : W^(-iV2^)0 = e"^ V ^ . (11.49) 
Definition II.4. The coherent vectors are the vectors of the form 

for some z & Z and the coherent states are the states of the form 

E.El . 

Definition II.5. A symplectomorphism T for the symplectic form Im(-, •) on a 
C-Hilbert space Z isa continuous R-linear automorphism on Z which preserves 
this symplectic form, i.e., 

\/zi, Z2 e Z : lm{Tzi,Tz2) = lm{zi, Z2) ■ 

A symplectomorphism T is implementable if there is a unitary operator Ut 
on 5+ (-2) such that 

e Z , UtW^(-z)U^ = W{Tz) . 
In this case Ut is a Bogolubov transformation corresponding to T. 

We recall a well-known parametrization, in the spirit of the polar decomposi- 
tion, of implementable symplectomorphisms. 

Proposition II.6. The set of implementable symplectomorphisms is the set of op- 
erators 



T — uexp[f] = u 
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where u is an isometry and f is an antilinear operator, self-adjoint in the sense that 
\/z,z' G Z, {z^fz') = {z',fz), and Hilbert-Schmidt in the sense that the positive 
operator is trace-class. Equivalently, there exist a Hilbert basis {'fj)jm of Z 
and {fij)ij G £^(N^; C) such that 

oo 

f = ^ (-, ifj) ifi and Mi, J G : fij = fj-^ . 

Proof. On the one hand, every operator of the form T = u exp[f] with u a unitary 
operator and r a self-adjoint antilinear operator is a symplectomorphism. Since a 
unitary operator is a symplectomorphism, and the set of symplectomorphisms is 
a group for the composition, it is enough to prove that exp[f] is a symplectomor- 
phism. It is indeed the case since, for all z, z' in Z, 

Im(e''2, z) = lm{e^'z, cosh(f)z') + lm.{e^z, sinh(f)z') 
= Im(cosh(f)e'"2;, z') + lm.{z' , sinh(f)e^z) 
= Im(cosh(f)e'"2;, z') — Im(sinh(f)e^z, z') 
= Im(e~'"e''z, z') . 

The implementability condition is then satisfied if we suppose f to be Hilbert- 
Schmidt. 

On the other hand, to get exactly this formulation we give the step to go from 
the result given in Appendix A in [7] to the decomposition in Proposition III. 61 
In Q an implementable symplectomorphism is decomposed as 

T = ue"^ , (11.50) 

where n is a unitary operator, c is a conjugation and f is a Hilbert-Schmidt, self- 
adjoint, non-negative operator commuting with c. It is then enough to set f = of 
to get the expected decomposition. To check the self-adjointness of f, observe 
that, for all z, z' in Z, 

{z',fz) = {z',rcz) = {rz',cz) = {z,crz') = {z,rz') . 

For the convenience of the reader we recall the main steps to obtain the de- 
composition in Eq. (III. 501) . First decompose T in its C-linear and antilinear parts, 
T = L + A, then write the polar decomposition L = u\L\. It is then enough 
to prove that \L\ + u*A is of the form e'^'". From certain properties of symplec- 
tomorphisms (also recalled in [7 |) it follows that the antilinear operator u*A is 
selfadjoint and |Lp + Iz = {u*AY. A decomposition of the positive trace 
class operator {u*A)'^ = J2j "^^^i^] '^^^^ orthonormal basis of Z yields 
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I-^I = + Using that — )■ one can study the operator \L\ 

and u*A on the finite dimensional subspaces ker(|L| — fxlz) which are invariant 
under u*A. It is then enough to prove that for a C-antilinear self-adjoint oper- 
ator / such that //* = on a finite dimensional space, there is an orthonor- 
mal basis {^pk}k such that f{(fk) = ^^k- The conjugation is then defined such 
that c{Yj f^k^k) = S Pkfk and f = sinli^^(Aj)l on that subspace. □ 



III Pure Quasifree States 



III.l From Quasifree States to Pure Quasifree States 

Let Z be the C-Hilbert space L'^{Sa,A x Z2). We make use of the following 
characterization of quasifree density matrices. 

Lemma III.l. The set of quasifree density matrices and pure quasifree density 
matrices, respectively, of finite photon number expectation value can be charac- 
terized by 

f E Z, I] a Bogolubov transformation, 

CeC\Z),C>0, \\C\\Biz)<l} 



pQd = ^Tif]^^w{-iV2f)WQn*vw{-iV2fy 

f E Z, V a Bogolubov transformation^ 

Proof. We only sketch the argument, details can be found in flUHl. It is not diffi- 
cuk to see that any density matrix of the form W{-iV2f)W ^^^^iJW{-iV2f)* 
is indeed quasifree. Conversely, if p G is a quasifree density matrix then it 
is fully characterized by its one-point function fp E Z and two-point functions 
ilp,0!p)- Moreover, W{—i\/2fp)* pW{—iy/2fp) E cOJ is a centered quasifree 
density matrix with the same one-particle density matrix, that is, the density ma- 
trix Wi~tV2fp)* pWi-tV2fp) corresponds to (0,7p - fpf;,&p - fpf;). Ob- 
viously, 7p — fpf* is again trace-class and dtp — fpf* is Hilbert-Schmidt. Now, 
we use that there exists a Bogolubov transformation U which eliminates cup, i.e., 
ll*W{V2fp/i)* p W{V2fp/i)lS corresponds to (0, 0). While this is the only 
nontrivial step of the proof, we note that if U is characterized by u and v as in 
Lemma HlOl then there is an involved, but explicit formula that determines u and 
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V. Again 7^ is trace-class because the photon number operator Nf transforms 

1 /2 

under U* to itself plus lower order corrections, W Nf\} = Nf + 0{Nj + 1). 
Finally, it is easy to see that (0,7p,0) corresponds to the quasifree density ma- 
trix r(Cp)/Tr[r(Cp)] with Cp := %(1 + ipY^. Following these steps we finally 
obtain 

as asserted. The additional characterization of pure quasifree density matrices is 
obvious. □ 

Lemma III.2. Let U G B{^) be a unitary operator. The following statements are 
equivalent: 

U G z^' a Bogolubov transformation; (III.51) 

<^3T implementable symplectomorphism, (III. 52) 

U = Ut, UTW{f)U* = W{Tf). 

G B{Z), V G C\Z) 'if eZ: (III.53) 

\]a*UW = a'iuf) + a{jvj f)- 

<;=4>1LJ = exp{iH), where H = H* is a semibounded operator, (III.54) 
quadratic in a* and a and without linear term. 

Proof. Again, we only sketch the argument. First note that (IIII.51I) <^ (IIII.52I) is 
the definition of a Bogolubov transformation. Secondly, UtW {f)U^ = W{Tf) is 
equivalent to f7r<l>(/)f/^ = $(T/). Hence, using that a* (/) = ^[$(/) -z<l>(i/)] 

and a(/) = ^[*(/) + i^{if)] we obtain the equivalence (IIII.52I) <^ (IIII.53I) . 
Thirdly, setting Ux = exp{i\H) and ■= Uxa*{f)Ul, we observe that 

dxa\{f) = i[H, a\{f)]. Furthermore, [H,a\{f)] is linear in a* and a if, and 
only if, H is quadratic in a* and a. Solving this linear differential equation, we 
finally obtain (IIlL53l) ^ (lnL54l) . □ 

Lemma III.3. For all Bogolubov transformation U and g G Z: 

W{g)U£ldV*W{g)* = OJ, (III.55) 
Uc0^?U* =c£ld- (III.56) 

Remark III.4. A pure quasifree state is a particular case of quasifree state with 
C = 0,thatisr(C) = m*. 

We come to the main result of this section. 
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Theorem IIL5. Let 0<a<A<oo, (^gM and p G M^, \p\ < 1. Minimizing 
the energy over quasifree states is the same as minimizing the energy over pure 
quasifree states, i.e., 

EBHF{.g,P,(rA) ■■= inf Tr[ifg,p-p] = inf Tr[ifg,p-p] . 

peas p6pa5 

For the proof of Theorem IIII. 5 1 we derive a couple of preparatory lemmata. 
Proposition III.6. Let C a non-negative operator on Z, then 

|Tr[r(C)] < ooj ^ ^CeC\Z) and \\C\\b[z)<^Y 

In this case Tr[r(C)] = det(l — C)^^ (We refrain from defining the determinant.) 
For the direction <^= the non-negativity assumption is not necessary. 

Proof. Let us decompose Z = 0j>oCej where C = with {ej)j>o an 

orthonormal basis of Z. Then d+{Z) = <S>j>o 'S+i^'^j) ^^d 

Tr[r(C)] = Tr[(g) T{c,)] = J] Tr[r(c,)] = J] 

j>0 j>0 j>0 

and the infinite product converges exactly when C G C^iZ) and ||C||/3(2) < 
1. □ 

Lemma III.7. Suppose Z^ is of dimension d < oo. Then, for any non-negative 
operator Cd such that Cd G C^{Zd) and \\Cd\\i3(Zi) < 1. there exist a 
non-negative measure fid (depending on C) of mass one on Zd and a family 
{pd{zd)} za£Za ofpure quasifrec states such that 



Tr[r(C)] 



Pd{zd) dfid{zd) ■ 



Proof. In finite dimension d we can use a resolution of the identity with coherent 
states (see, e.g., yj) 



lr(2,) = \ Ez,K, 



where Zd is identified with and dzd = dxd dyd, Zd = Xd + iyd- Using Equa- 
tion (III.49I ) we get 

dzd 



T{C) = I T{C'/')E,,EIT{C''' 

eM\C'"^d?-\zd?)dzd 
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The measure dudizd) = tt '^expdC^/^Zrfp — \zd\^)dzd/T'c\J^{C)] has mass one. 
Indeed 



f exp(-z^(l^^ - C)zd)^ = f\[ 



2 2\\d'^d,y 



{-{1 - c,){x' + y')) 



TT 



nr^ = Tr[r(C)] 



where C = J2'j=i "^^^^ (^i)i=i orthonormal basis of Zd- □ 

Proof of Theorem 51 The inclusion pllS^ C 05^ implies that 

inf Tr[ifg,p^p] < inf Tr[i7g,p^p], 
peas pepQA 

and it is hence enough to prove for any quasifree state 



that the inequality 



TilHg^j^Pgf] > inf Tr[if3,p-p]. 



holds true. The operator C is decomposed as C = '^j^QCjeje* where (cj) is an 
orthonormal basis of the Hilbert space Z and cj > 0. Let Cd = J2j<d Let 

then using Lemma |IIL7] with Zd = ®j<d Ccj, ^+2 = d+{Zd © Z^) = d+Zd ® 
d+2d and the extension of the operator r(Cd) on d+2d to ^^Zd^d+Zj- by T{Cd)® 
(il^^ifi* j^) (which we still denote by T(Cd)), we obtain 

Pqf,d = / Pd{Zd) dpd{Zd) , 



d 



where pd{zd) are pure quasifree states and the pd are non-negative measures with 
mass one. Note that 

' ■ Tr[r(c)] n;' ^ 
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as d — > oo. Further note that p^/ > Vd Pqf,d, for any d e N, since V{C) > V{Cd). 
Thus 

^ Vd j TY[Hg^ppd{zd)\ djidizd) 
> Ud inf Tr[Hg^j;pd{zd)] 

for all d gN, and in the limit d — > oo, we obtain 

Tr[Hg,pP,f] > jini{^.} mi Tr[Hg,,^p] = inf IV[i/,,^-p] . □ 

III.2 Pure Quasifree States and their One-Particle Density Ma- 
trices 

Let Z be a C-Hilbert space. 

Definition III.8. Let p E DDJl be a density matrix on the bosonic Fock space 
over Z. If Tr[pN/^] < oo, we define p^-^ e through 

V(/?, V e -2 , V*''V'V' = Tr[a*{ipya{jPYp] . 

We single out 

i.e., fp E Z is the unique vector such that Tr[a{ip) p] = ip* fp, for all ip E Z. 
Furthermore, with p — W{y/2fp/i)*pW{y/2fp/i), the matrix elements of the 
(generalized) one-particle density matrix are defined by 

7p = p''' e B''' and ap = p^'^ e S^'^ ^ , 

in other words 

V(^,'0eZ: (^,7p^) = Tr[pa*(v9)a(^)] , 
® , ttp) = Tr[pa(^/;)a(<^)] . 

Note that fp, 7p, and ap exist for any p e since A^/p, pNf e 

Remark III.9. For a centered pure quasifree state p, p^'^ vanishes when p + g is 
odd. 
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Remark III.IO. Another definition of the one-particle density matrix 7^ would be 
through the relation 7p(y9) = Tr [a* {(f) a{ip)p]. We prefer here a definition with 
a "centered" version p of the state p, because this centered quasifree state p then 
satisfies the usual Wick theorem. The same considerations hold for ap. 

Hence, any quasifree density matrix is characterized by (/p, 7p, ap), since p^''^ 
can be expressed in terms of (/p, 7p, ap). 

When fp = 0, the definition of 7p is consistent with the usual one, for zi, 
Z2 e Z, {zi,'^pZ2) = Tr[a*(^2)a(^i)p]- The definition of ap is related with the 
definition of the operator dp (here denoted with a hat for clarity) used in the article 
ofBach, LiebandSolovej [|3J, through the relation (zi®Z2, ap).E'»2 = {zi,apCZ2)z 
with c a conjugation on Z. 

Example III.ll. A centered pure quasifree state satisfies the relation, 

p2:2 = ^(g)^ + ^(g)^Ex +aa* G -B^'^ , (III.57) 
where the exchange operator is the linear operator on Z®^ such that 

'izi, Z2 G Z, Ex(2;i (g) Z2) = Z2® zi 

and where for any h G Z®'^, aa*b = {a, b)z<s2 a. 

We now tum to another parametrization of quasifree states, by vectors in a real 
Hilbert space. This parametrization enables us to use convexity arguments. 

Proposition III.12. Let T = ue^ be an implementable symplectomorphism and p 
a quasifree state of the form p = 1]^QQ*Vt- Then 

-fp = |(cosh(2f) - 1) , (III.58) 

\/zi,Z2 e Z : {zi® Z2,ap)z^2 = {zi,lsmh.{2f)z2) . (III.59) 

Proof of Proposition 1 21 We have Ti = ue^i = uie~^ = iue~^ and for all 

z e z 

TT[pW{-iV2z)] = TT[l]*Tnn*UTW{-iV2z)] 
= Q*W{u/{-iV2z))Q 
= n*W{-iV2ue-^z)n 
= exp ( — ||ue^^2;|^) 
= exp ( — i|e"''2;|^) 

From this formula we can easily compute the function 

h{t,s) := TT[pW{-UV2z)W{-siV2z)] = exp ( - ||e-^(t + s)^!^) 
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whose derivative dtdg at {t, s) = (0,0) involves a and 7 

dtdsh{0,0) =Ti[p{a*iz) - a{z))^] 

= -2z*-fz + 2Re(a*2^2) - z*z . 

But we have also 



dtdsexp{-l-\e '~(t + s)z\'^] 



2' 



t=s=0 



(e-'zYie-^'z) 



-{cosh{f)z — sinh(f)2;)*(cosh(f)z — sinh(f)z) 
-(cosh(f)2;)*(cosh(f)z) 
+ 2Re(sinh(f)2;)*(cosh(f)2;) — (sinh(f)2;)*(sinh(f)2;) 
-2;*(cosh^ f + sinh^ r)z + 2Re(z*(sinh f cosh r)z) 

-z* cosh(2f)z + 2Re(z*- sinh(2f)2;) 



and hence, using the polarization identity 



4zVz' = (z + zT' - [z 



Z 



to recover every vector from from linear combinations of vectors of the form 
z"^^. we arrive at (inL58])- (IIIL59l) . □ 

Proposition III.13. The admissible 7, a for a pure quasifree state are exactly 
those satisfying the relation 

7 + 7^ = (" ® 1)* (1 ® ") , (ni.60) 

with 7 > 0. 

This is the constraint when we minimize the energy as a function of (/, 7, a) 
with the method of Lagrange multipliers in Section IVIIII 

Proof. If 7, a are associated with a quasifree state, then there is an f such that 7, 
a and f satisfy Equations (IIII.58I) and (IIIL59I) . then 

{zi, {a* ® 1)(1 ® a)z2) = {a* ® zl){z2 ® a) 

= {[a*{z2 ® 1)] ®zl)a 
= {a*{z2 (S) l),isinh(2f)zi)2 
= (a*, 2:2 ® I sinh(2f)zi)^^2 
= (isinh2(2f)zi,Z2>^ 

= ((i(cosh(2f) - 1) + i(cosh(2f) - 1)2)^1, Z2)^ . 
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Conversely, if 7 and a satisfy Eq. (1111.601) then we define the C-antilinear operator 
a such that (zi, az'^ = {zi (g) Z2)*a, and set f = ^ sinh~"'^(2d), then 

Wzi,Z2eZ: {zi® Z2,ap)z2 = {zx^az^) = (zi, | sinh(2f)z2) , 

which, in turn, implies that {a* (g) 1)(1 a) = ^ sinh^(2f). Hence, we have 

7 + 7^ = - sinh^(2f) 

and as 7 > 0, it follows that 7 = |(cosh(2f) — 1). Then 7, a is associated with 
the centered pure quasifree state whose symplectic transformation is exp[f]. 

□ 

IV Energy Functional 

— * — * 

Notation: We first recall that, as before, we denote by k, and \ k\ the multiplication 
operators k 1^2 and \k\ lc2 on Z = L?'{S„j^ x Z2), with three components in 
the case of k. 

We now work at fixed values of total momentum p eM?. The operator Hg^^ is 
given by 

Hg,p = l{dT{k) + 2Re a* (G) - p)'^ + dT {\k\) , 

where G{k) = G{k, ±) := g\k\~^/'^e±{k). The energy of a pure quasifree state p 
associated with f E Z,^ ^ ^^{Z), a G is 

^5,p(/,7,«) := "^AHg^pP], (IV.61) 

where Z is the C-Hilbert space Z = L'^{S^^a x ^2) and C^{Z) is the space of 
trace class operators on Z. 

Proposition IV.l. The energy functional (IIV.61I) is 

£,,p^(/,7,«) = i{(Tr[7^] + r^/ + 2Re(rG) - p}'' 
+ Tr[7fc-7fc] +a*{k-^k)a + TT[\k\^^] 

+ 2Re{a*[(G + kf)-'"']} + Tr[(27 + + kf) ■ {G + kf)*]} 
+ Tr[7|fc|] + r|fc|/. (IV.62) 
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where the following positivity properties hold 

(Tr[7fc] + rkf + 2Re(r G) - p)'^ > , 
Tr[7A^ ■ jk] + Tr[7fc]-2 + a*{k- ®k)a + Tr[|A?p7] > , 
(Tr[7^] + rfc/ + 2Re(rG)-p)-2 
+Tr[7A; ■ 7^] +Q*{k- ®k)a + Tr[|A;p7] > , 
2Re(a*((G + fc/)'^')) + Tr[(27 + 1)(G + fc/) ■ (G + kf)*] > . 

The energy of a pure quasifree state in the variables f and r is 

kAf^ r) = H(Tr[|(cosh(2f) - 1)^] + f*kf + 2Re(/*G) - p)'' 
+ Tr[i(cosh(2f) - ■ i(cosh(2f) - l)k] 
+ Tr[i smh{2f)k ■ | sinh(2r)fc] + Tr[|fcpi(cosh(2f) - 1)] 
+ 2Re(i sinh(2f)(G + kf); {G + kf)) 
+ Tr[(2|(cosh(2f) - 1) + 1)(G + kf) ■ (G + kf)*]} 
+ Tr[i(cosh(2f) - 1)|A^|] + f*\k\f . (IV.63) 

Proof. Using the Weyl operators, 

£g,p{f^^^(^) ■= ^AHg,pp] = Tr[/fg,p.(/)p] 

where H^^) = W{V2f /iyHg^^W{V2f /i) and p = W{V2f /if pW{V2f /i), 
so that p is centered. Modulo terms of odd order, which vanish when we take the 
trace against a centered quasifree state, Hg^p{f) equals 

HgAf) = lid^ik) + f*kf + 2Re(a*(fc/ + G)) + 2Re(/*G) - p)'' 
+ dT{\k\) + f*\k\f + odd 
= ^{dVik) + f*kf + 2Re(rG) -p)-' 

+ ^{2Re{a*ikf + G))}' + dT{\k\) + f*\k\f + odd . 

To compute £{f, 7, a) we are thus lead to compute, for E and u G M'^, 

Tr[p(dr(fc) and Tr[p (2Re{a((^)})-2] . 

The expression of the energy as a function of (/, 7, a) then follows from Propo- 
sitions |IV2] and |lV4l The expression of the energy as a function of (/, r) follows 
from Proposition UlL 1 2[ □ 
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Proposition TV.2. Let u G M^, then 

< Ti\~p{dT(k) + = (Tr[7A^] + n)-^ - Ti[-ik]-^ 

+ Tr[7A^ ■ 7A;] + Tr[7A;]-2 + a*{k-0k)a + Tr[|A^p7] . 

This condition is used with u = p- f*kf - 2Re(/*G). 
Proo/ Indeed, 

(dr(fc) + u)-2 = rfr(fc)-2 + 2cir(A?) -M + u'^ 

Then we use that Tr[pdT(k)] = Tr[7A;], add and substract Tr[7A;]'^ to complete 
the square and compute Tr[p dT{ky'^] using Lemma lTV.3[ □ 

Lemma TV.3. Let X G B^'^, then 

< 1i[pdV{X)dV{Xy] = Tr[7X7X*] + |Tr[7X]|2+a*(X(g)X*)a+Tr[XX*7] . 

Proof. Indeed, using Equation (IIII.57I) . 

TT[pdT{X)dT{Xy] 

= Ti[p{Jx{ki, k[)X{k2, k'2)a*{ki)a*{k2)a{k'2)a{k[)dkidk2dk[dk'2 + dT{XX*)] 

= Tr[(7®7 + 7(g)7Ex + aa*)(X®X*)] + Tr[7XX*] 

= Tr[7X]Tr[7X*] + Tr[7X7X*] + a*{X (g) X*)a + Tr[7XX*] . □ 

Proposition IV.4. Let G Z, then 

< TT[p{a*{ip) + a{ip)y] = 2Re(a*(v;^2)) ^ rp^^^ ^ ■^^^^*^ (jy_g4) 

and \2Re{a*{if'^^))\ < Tr[(27 + l)ipip*]. 

This condition is used with the three components of ip = G + kf . 

Proof. A computation using the canonical commutation relations yields 

Tr[p{a\^) + a{^)f] 

= Tr[p {a*{^)f + p {a{p,)f + p {a*{p)a{^) + a{p,)a*{^))] 
= + <p'^^*a + Tr[7(^(^* + (7 + 1)^^*]- 

□ 
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V Minimization over Coherent States 



For this section we can take a = if we consider the parameter / in the energy 

to be in Z := L'^{S^^a x Z2, + \k\)dk). Recall that S^^a = {k eR^\a < 

\k\ < A}. 

Remark V.l. For a coherent state (see Definition 111.41) the energy reduces to 

^aAf) = ^ll^ir + lifkf + 2Re(rG) - p}'' + /* + \k\)f . (V.65) 

Note that, for a > 0, Z = L'^{Sa,A x ^2, dk) = Z, while for cr = 0, Z C Z, and 
£g,pif) extends to Z by using Equation (IV.65I) . 

Theorem V.2. There exists a universal constant C < 00 such that, for < a < 
K < 00, ln(A + 2) < C and |p| < 1/3, there exists a unique which mini- 
mizes £g^^ in Z. 

1. The minimizer solves the system of equations 

fp = , ^'i^ ^ ^ , (V.66) 

u^ = V- 2Re(/;G) - , (V.67) 

with \u^\ < \f>\. 

2. ForO <a < A< 00, 



jnf ^3,p(/) = inf ^<;,p(/) = £gAf} 



pj } 



fez 

and for < a < A < 00, we have that G Z. 

3. For fixed g, a, A, as a function off), 

4. For all f in Z, 

SgAfp + f) = SgMp) + nilM' + \k\ - ur^)f 

+ \{f*kf + 2Re{flkf) + 2Re(/*G))-' . (V.68) 
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5. The energy Sg^f{f^) of the minimizer compared to the energy of the vacuum 
state Sg^^{0) is 

£gAfp) = ^:,,p(0)-i2Re(/;%-G)-i|Mp^-pl2_ 

Note that the term 2Re{ f^u^ ■ G) is non-negative. 

Remark V.3. Our hypotheses are similar those of Chen, Frohlich, and Pizzo (HI, 
where their vector V EZ is analogous to up'm our notations. 

The construction of as the solution of a fixed point problem and the depen- 
dency in the parameter p imply that the map p^u^h of class C°°. 

Remark V.4. We note that we also expect to have up in the neighboorhood of p. 

Remark V.S. The minimizer is constructed as the solution of a fixed point prob- 
lem. As a result the application 

(c^, A, ^, p) ^ inf TilHg^pp] 
pecoh 

is continuous on the domain defined by Theorem IV. 21 and at a, A fixed, 

{g,p)^ inf TT[Hg^pp] 
pecoh 

is analytic for g"^ < C/ ln(A + 2) and \p\ < 1/3. 

Proof of Theorem \V2\ Proof ojU\ Assume there is a point fp where the minimum 
is attained. The partial derivative of the energy at the point fp 

= {{ffkfp-p + 2Re(/|G)) ■k + ^\k\' + \k\)fp- {p- flkfp- 2Re(/|G)) • G 

then vanishes, where the derivative df*£{f) at a point / is the unique vector in 

Z* = L\S^,A, + \k\)-^dk) defined by 

£{f + 5f) - 8U) = 2Re(r df,S{f)) + o(||5/||^) 
with /, 5f E Z. Observe that 

< £g,piO) - ^gAfp) 
= hp\-' - ^(Wp + 2Re(/;G) - p)-' - f;ihk\' + \k\)fp 
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and hence \p\ > \u^\ with := p — f^kf^— 2Re(/iG). Since \up\ < |p| < 1, it 
makes sense to write 

- up-G 

"21^1^ ~\~ I k I ' k 

Hence the minimum point satisfies Equations (IV.66I) and (IV. 67 1) . It is in partic- 
ular sufficient to prove that there exist a unique up in a ball 5(0, r) with r > |p| 
such that the function in Equation (|V.66I) satisfies also Equation (|V.67I) to prove 
the existence and uniqueness of a minimizer. 

Proof of the existence and uniqueness of a solution. Let |<r<l,MGM^, 
|u| < r < 1 and 

^\k\^ + \k\ - k ■ u 

Observe that G Z, indeed, if \u\ < 1 then ||A;p + -k-u> (1 -r)(i|fcp + 
\k\), and with e{k) = e{k, +) + e{k, —), 



/ 

J\k 



:l\k\' + \k\)\Mk)\'dk <g^ I ^ ^ ^^i^rif^in ^fe < +^ • 



<Co/ln(A + 2) 



2 1 
2 



(1 -r)2 

for some universal constant Co > 0. Observe then that 

r \Gik)\'^ dk<Cog'HA + 2) 

J\k\e[a,A] II^P + 1^1 

for some universal constant Co > 0. It follows that $tG G L^{Sa,A x Z2). Note 
that if a = then <l>^ ^ ^^(S^^ x Z2) (for m ^ 0). 
We can thus define the application 

5(0, r)3u^ ^(m) ■.= p- - 2Re($tG) G . 

We check that the hypotheses of the Banach-Picard fixed point theorem are veri- 
fied on the ball B(0, r), which will prove the result. 

Stability: If g'^ ln(A + 2) is sufficiently small, we get from 

\^iu)\ < \%k<!>^\ + \2Rei%G)\ + \p\ 

and the estimates above that the sum of the two first terms is smaller than r — 1/3 
and since \p\ < 1/3 the map ^ sends 5(0, r) into itself, 

^(5(0, r)) C 5(0, r). 
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Contraction: For u and v in B(0, r), we have that 

\^.{k) - <^,^{k)\{l\k\' + \k\) 

:\m' + \k\) 



u.G{k) v.G{k) ,,1-^2 



\\k\'^ + \k\-k-u \\k\'^+\k\-k.v 

u-v\\G{k)\ 



2|A;|2+ |A;| -k-u 

1 



+ \v\\G{k)\ 1-^ ^ ^ ^ ^1 + |fc|) 

< 1-^(1 + — — ^ffl ) 

(1-0 \\k\'+\k\{l-ry 
<\u-v\\G{k)[ ^ 



(l-r)2 ■ 

For the term 2Re($^G), we observe that 

|2Re($tG)-2Re($tG)| 

< Ci/ln(2 + A)2|ir-'y| ^ 



(l-r)2 ■ 

Note that, for g"^ ln(2 + A) < (1 - r) V(3Ci), 

|2Re($:-G)-2Re($:-G)|<i|K-i;| 

Finally, for the term $lA;$jj, we obtain the estimate 

\%k^a-^%k^^\ 



< f dl^l' + \k\)\Mk) - MkMMk)] + \Mk)\)d'k 

J\k\e\a,A\ 

[ \Gik)\i\^4k)\ + \%^{k)\)dk 

J\ma,A] 



'|fe|ek,A] 
\u — v\ 



(1 -r)2 



\k\e[<j,A] 

hu 

(1 -r 



2 1 



< C2\u-v\ {\u\ + \v\)g'H^ + 2) , 
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and thus this term can be controlled for |(7 ln(A+2) p sufficiently small by | |m— | . 
We thus get a contraction 

2 

- ^{u)\ < -\u-u\ 

and with = Equation (IV.66I) is solved. 

Proof of \3\ The expression of the energy £g,p{f) given in Equation (IV.65I) 
implies that £g,p{f) > \ HGp, and for p = this minimum is only attained at the 
point = 0. It follows that f.^ = d^f^ ■ p + . From Equation (IV. 67 1) we 

deduce 



u 



p 



p-2Re{idp^fg-p}*G) + Oi\p\ 



2\ 



and thus 



_ ip-2Re{idp^f^,.p)*G)).G , 
l\k\2 + \k\ - k ■ Uj, ^'^'^ 

= + \k\)-\p- 2Re{{dp^f^ ■ p}*G)) ■ G + 0{\p\') . 

Expanding the left hand side of this equality in brings 

dpfo ■ P= il\k\' + \k\)-\p- 2Re((9,^/5 ■ p)*G)) ■ G 

and hence d^f^^ = {^\k\'^ + \k\ + 2G ■ G*Y^G. The expansion of f^io the second 
order is then 

U = il\k\' + \k\ + 2G ■ GT'G ■ p + 0{\p\') . 

We can compute the energy modulo error terms in To have less heavy 

computations we set A = ^|A;p + \k\ + 2G ■ G* and get 



2" " 2' 

^ bf + ^(2Re(p ■ dj,fiG)-p)-'+p- dj^fiihkl' + \k\)dp^fi.p 



= -(2Re(p ■ G*A-^G)y^ - 2p- G*A-^G ■ p 

+ p.G*A-\]^\k\' + \k\)A-^G-p 
= 2{p- G*A-^Gy^+p- GM-^((^|A?p + \k\) - 2A)A-^G ■ p 
= p- G*A-^2G ■ G*A-^G ■ p- p- G*A-\]^\k\^ + \k\ + AG ■ G*))A-^G ■ p 
= -p- G*{\\k\^ + \k\ + 2G ■ G*)-^G ■ p 
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which yields the result. 

Proof oj^ The Taylor expansion of the energy around /p-is 

+ l{{rkf + 2Re(rG) + 2Re(/|fc/))-2 

+ 2(/;%+ 2Re(/;G) - p) ■ fkf + r\k\'f} + r\k\f . 

Since df*S{f^) vanishes this gives Equation (IV.68I) . 

Proof of\^ It is sufficient to replace / by — /p- in Equation (IV.68I) . The obser- 
vation 

2" I ^ I ^ ~l~ I ^ I ^ ' 

shows that 2Re{ f^up ■ G) is non-negative since l^p-l < 1. □ 

VI The Minimizer for the Energy Functional vary- 
ing over Pure Quasifree States 

Definition VI.l. Let Z be a C-Hilbert space. Let Y be the M-Hilbert space of 
antilinear operators f on Z, self-adjoint in the sense that Vz, z' E Z, {z, rz') = 
{z',fz), and Hilbert-Schmidt in the sense that the positive operator is trace 
class. The space X = Z x Y with the scalar product 

{{f,f),{f,f'))x = rf + TT[ff'] 

is an M-Hilbert space. 

Keeping a > 0, we only need to use Z = L'^{Sa,A x Z2) in this section. 

Theorem VI.l. Let < a < A < 00. There exists C > such that for g, \p\ < C 
there exists a unique minimizer for £g,p{f, f). 

Proof. This result follows from convexity and coercivity arguments. By Proposi- 
tion lVL3[ £g,p{f, f) is strictly 6'-convex (i.e., uniformly strictly convex) on Bx{0, R) 
for some R > and 9 > 0. Since £g,p{f, r) is strongly continous on the closed 
and convex set -Bx(0, R) of the Hilbert space X we get the existence and unique- 
ness of a minimizer in Bx{0, R). (See for example [HI. The uniform strict con- 
vexity allows to prove directly that a minimizing sequence is a Cauchy sequence.) 
Proposition lVL4l then proves that it is the only minimum of £g,p{f, f) on the whole 
space. 

Note that to use Propositions IVL3I and IVL4I we need to restrict to values of g 
and |p| smaller than some constant C > 0. □ 
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Proposition VI.3 (Convexity). There exist < C, R < oo such that for g < C 
and |p| < |, the Hessian of the energy satisfies l-LSg^^f,r) > ^Ix on the ball 

Proof. We use that strict positivity of the Hessian implies strict convexity and thus 
first compute the Hessian in (0, 0). The Hessian HSg^p{f, f) e B{X) is defined 
using the Frechet derivative 

= DSgM,r){Sf,Sf) + l{{Sf,6f) , nSgM^r) {Sf,r))^ + o(||(5/, 5f)|||) 

with DSg^^O, 0) e B(X,R). (Note that differentiabiUty is granted in this case 
because |^| < A < oo.) For any /i > 0, V(/, r) e X, 

((/,f), ^^4p(0,0) {f,f))x 

= 2Re{rkf; G) + ^(2Re(/*G))-^ + Trff^G • G*\ 

+ ^{Tr[rA; -fA-] + Tr[|A;|V]} 

+ Tr[f2^|^| -k-v)\-r r (^I^P -r\k\-k-v)f 

A* 

+ ^{(2Re(5/*G))-' + Tr[ffc-ffc] +Tr[|fc|2f2]| 

+ iv[f2(i^i - ^ • p)] + r (^i^r + 1^1 - ^ • v)j 

> Tr[f\\k\ -k-p+{l- ii)G ■ G*)] + r ((i - -)\k\' +\k\-k- p)f , 

Zi jjj 

since 

|2Re(f^/;G)| < 2||fG|| ||^/|| = 2V;[I||fG||^||^/|| < /.||fG||2 + -Wkff . 

V A* A* 

With = 2 we obtain (with |p] < i) 

^7{4p-(0, 0)(/, f) > TV[f^(|^| -k-p-G- G*)] + f*{\k\ -k-p)f 

> Tr[f2(|^|(l - \\\k\-'^'Gr) -k.p)] + f*{\k\-k. p)f 

>iv[fV(l-iii^i-v2G||^)]+r|/ 
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and for g small enough 

Ins, AO, 0) > ^ . 

We then compare it with the Hessian in points near zero. Observing that the 
Hessian is continous with respect to (/, f, p, g), we deduce that there exist R < oo 
and C > 0, as asserted. □ 

Proposition VI.4 (Coercivity). Suppose p and C > are fixed such that \ \p\^ -\- 
< aR^, with the value of R given by Pwposition \VI.3\ for any < g < C. 
For every {f,f) G X, 

4p(/,r) > Tr[f2|^|] + r|^|/ > cr\\{f,r)\\l. 

Since SgA^^O) = + < o-R^, any minimizing sequence takes its values 

inBx{0,R). 



VII Asymptotics for small Coupling and Momen- 
tum 

We use below an identification between self-adjoint C-antilinear Hilbert-Schmidt 
operator f and symmetric two vector r given by the relation (yj, rip)z = ® 
if), r)zm. Note that the self-adjointness condition for f is equivalent to the sym- 
metry condition r G 2^^. 

Theorem VII.l. Let < a < A < oo. There exists C > such that for 
\g\, \p\ < C, there exist two functions f g^p and fg^p which are smooth in {g,p) such 
that the minimum of the energy £gAf-> ^) '■^ attained at {fg^p, Tg.p)- These functions 
satisfy 

Up= ,r^ +o{Ug,p)r) 

rg^.= -S-'G-WG + 0{\\{g,p)f), 
with S = k ■ ®k + 2(||A;p + \k\) V Iz- As a consequence 

EBHF{.g,P,cr,k) 

= SgAOx)-p- G \ ^ G-p- \g-^'*S-'G-^' + 0{\\{gM') ■ 
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Remark VII.2. The energy in Ox is the energy of the vacuum state and is £g,p{Ox) 



+ IIIGp. Further note that 



\k\^ + \k\ 



(G-p) = g'\p\'[2n'-- In 



in 



A + 2 
a + 2 



and in particular does not depend on the choice of the polarization vectors e. 

The quantity G^^* 5*^^(5^^ does not depend on the choice of the vectors e 
either since 



,.u=±-' ^\k,\\k2\S{k,M) 

and with is the orthogonal projection on u in M?, 



1 + 



ki 

ki k2 



Proof of Theorem WII. 1\ Let 
and(/)(^7,p):= ( ) such that 



then a derivation of Equation (IVII.69I ) with respect to (/, r) brings 

The term which is independent of {g,p) and quadratic in ( ^ ) in the energy is 



(VII.69) 



-{Tr[f5f] + r(|fcp + 2|fc|)/} 



thus, in (Og,p-, 0/,r), 



df,rF 



\k\^ + 2\k\ 

s 
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To compute dg^j;F in 0, observe that no part in the energy is linear in {g,p) and 
linear in (/, f). Thus dg.i;F{0g.^,0f.^) = and we get 

5<;,p/(0,,p-)=0. 

Differentiating a second time Equation (IVIL69I ) brings 

= dl,F + 2df,,dg,,F O dg,, ( / ) + dj,,F O ( / ) + dlF{dg,, ( ^ ) , 9, ( ^ ) ) . 

Since dg^^ ( I ) (Og^p) = 0, it follows that 

( i ) (Og,^) = -[df,rFiOg,p^, Of,r)r'dl^^FiOg,p^, 0;,,) . 

The part of the energy which is quadratic in {g, p) and linear in (/, f) is — 2Re(/*(j) • 
p + Re(fG; G), it follows that, in (Oc,,p-, 0/,^), 

which gives in Og^p- 

V-s-'la^G o)-v(a,G 0) 

Hence the expansion of ( ^ ) up to order 2. 

We can thus express the energy around O^^p- modulo error terms in (9( || || ^) 

niin4,p(/,r) - 4,p(0, 0) 

= \{{i:i\f% + rkf + 2Re(/*G) + Tr[ffc ■ ffc] + Tr[|A?pr2] 

+ 2Re(f(G + kf)- {G + ^/)) + ||Gf + /* |^| V} 
+ Tr[f2|^|] + r|fc|/-4p^(0,0) 

^ -2Re(rG) • p+ |Tr[f5f] + Re(fG; G) + /* (^l^p + \k\)f 
^ -2Re(rp • G) + |Tr[f^f] + Re(fG; G) + /* + |^|)/ 

^ ^ (P-G)*(P-G) ^ 1^.V2*^-1^.V2 _ ^.V2*^-1^.V2 _^ (P ■ G)*(P ■ G) 

2 \\k? + \k\ 

which completes the proof. □ 
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VIII Lagrange Equations 



This section formulates the results of Section |Vl]in terms of 7 and ol subject to the 
constraints 7 + 7^ = [a* ® lz){lz ® ct), without reference to the parametrization 
of 7 and a in terms of r. 

Suppose f e Z,a e Z'^\ 7 G C\Z), X e B{Z) = Bandu e Let 
^(A) = iA; ■ VA; + A V 1 and ^(7) =7 + 7 
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Theorem VIII.l. Suppose (/, 7, a) is a minimum of the energy functional £ such 
that ||7||b(^) < \. Then there is a unique (A, u) such that (/, 7, a, A, u) satisfies 
the following equations, equivalent to Lagrange equations 

M{-f,u)f = -(% + ^1) -m) -G-k-yiG + kfYa (Vm.VO) 
A{X)a = + A^/)-^2 (VIIL71) 

7 = G'\{a* ® ® a)) (VIIL72) 

POO 

A = / e-*(5+7) (M(7, u) + {G + kf) ■ {G + kfy)e''^'^+^^dt 
Jo 

(VIIL73) 

u=p~ Tr[7fc] -f*kf- 2Re{ f*G) (VIIL74) 

with M(7,'u) = + \k\ — k ■ u + k ■ ^yk. 

Assuming |p| < sufficient conditions such that M(7, u) and A{\) are in- 
vertible operators are |u| < 1/2, 7 > and || A — {\k\'^/2 + \k\ — p- k)\\is < a/2. 
Equations AVIII.70\) to Will. 74\ then form a system of coupled explicit equations. 

Remark VIII.l. To prove that Equations (IVIII.70D to (IVIIL74D admit a solution 
we use here the result of existence of a minimizer proved in Section|Vll It can also 
be proved directly by a fixed point argument by defining the applications 

*^(/, a, 7, u) = -M(7, u)-\k{^ + U.)-u)-G-k-y{G + kfYa 

vI/„(/,A) = -^(A)-ii(G + A^/)-^^ 
^j{a) = g-\{a* ^ lz){lz ^ a)) 

roo 

^a(/,7,^)= / e-*(^+^)(M(7,n) + (G + ^/)-(G + ^/)*)e-*(^+^)cit 
Jo 

^uif, 7) = P- Tr[7^] - f*kf - 2Re(r G) 
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defined on balls of centers centrers 0, 0, 0, ^\k\'^ + \k\ — k.p and p and proving 
that the application 

vl/(;,,)(/, \) = {^f [/, vI/„ {/, A} , VI/, {vI/,(/, A)} , ^^{f, ^,{^^ [/, A])}] , 

vl>, [/, {vl>„(/, A)} , vl>,^ {/, [f, A])}]) 

is a contraction for a convenient choice of the radiuses and a sufficiently small 
coupling constant g. Note that it is then convenient to consider the norm of 

L\S^,AxZ2,\k\^)forf. 

Proof of TheoremMUJ} Indeed, set u = p - Ti[-fk] - f*kf - 2Re(/*G) and 
define the partial derivatives as df*£{f,'y,a) E Z, 7,a) G Z'^'^ and 

d^£{f, 7, a) e B{Z) ^ C}{Z)' such that 

^(/ + 5/,7 + 57,« + H-^(/,7,«) 

= 2Re(r 7, «)) + 2Re(fc* 7, «)) 

+ Tr[57(9,£:(/,7,a)] + o(||(5/, ^7, "^a) Il2x£i(-Z)x2^2) . 

Recall the energy functional is given by Equation (IIV.62I) and this yields 

df*E{f, 7, a) = ^{2ikf + G) ■ (Tr[7^] + f*kf + 2Re(/*G) - p) 

+ 2k-V{G + kfYa + ■ (27 + 1)(G + kf)} + \k\f 
= -ikf + G)-u + k-y{G + kf)*a + ^ ■ (7 + + kf) + \k\f 

= M{-f, u)f + (% + ^1) - ■ G + ^ ■ V(G + , 

7, a) = . + ^(G + ^ 

d,EU, 7, «) = ^{2^ ■ (Tr[7^] + /*^/ + 2Re(r G) - p) 

+ 2k-^k+ \k\^ + 2{G + kf) ■ (G + kf)*} + \k\ 
= M{^,u) + {G + kf)-{G + kfr. 

The constraint given by Equation (IIIL60I) can be expressed as 

C(/,7,a) = (VnL75) 

with 

C : Z X C\Z) X ^ £1(2) 

(/, 7, a) ^ 7 + 7^ - (a* (g) ® a) . 
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Equation (IVIIL75I) is equivalent to Equation (IVIII.72I) . The application C has a 
differential DC{f, -f,a): Zx C\Z) x ^ ^i^^) ^^ch that 

DC(/,7,«)(5/,57,M 

= 57 + 57 7 + 7(57 - ((5a* l2)(l2 ® a) - (a* ® lz)iiz ® <5«) • 

For ||7||b(2) < I the application DC(/, 7,0;) is surjective. Indeed it is already 
surjective on {0} x C^{Z) x {0}, since, for every 7' G the equation 5'y + 

(57 7+7 (57 = 7' with unknown 5'j has at least one solution, see Proposition I VIII. 3 [ 
We can then apply the Lagrange multiplier rule (see for example the book of 
Zeidler IfTOl ) which tells us that there exists a A G B{Z) such that 

V((5/, 6a, (57) , D£{f, a, 7)((5/, 6a, 6j) + Tr[DC(/, a, j){6f, 6a, 6j) A] = , 

that is to say 

2Re{6rdf,£{f, 7, a) + 6a*d^.£U. 7, «)) + Tr[.9^£(/, 7, «)(57] 
+ Tr[(57 + (577 + 7(57- ((5a*®l2)(l2®a)- (a*® l^)(l^®(5a))A] = 0. 

This is equivalent to Equations (IVIII.70I) . (IVIII.71I) and 

K\ + 7) + + 7)A = M(7, u) + {G + kf) ■ {G + kfr (VnL76) 

Using again Proposition IVIIL3I we get that Equation (IVIII.76I) is equivalent to 
Equation (IVIIL73I) . 

For the invertibility of ^(A) note that 

^(A) = \ik^i + i^ kr + (1^1 - k.p + x-l\k\^- \k\ + k.p)yi 

>(|-A-(|^|V2+|^|-p-^)b)lVl. 

For M{'y,u), M{'y,u) = + \k\ - k ■ u + k ■ 'jk > o-/2 if 7 > and 

|m| < 1/2. □ 

Let us recall a well known expression for the solution of the Sylvester or Lya- 
punov equation. 

Proposition VIII.3. Let A and B be bounded self-adjoint operators on a Hilbert 
space. Suppose A> al with a > 0. Then the equation 

AX + XA = B 

for X a bounded operator has a unique solution Xa{B) = e~^^Be~^^dt. 
If B a trace class operator then the solution X is also trace class. 
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Proof. Indeed, Xa{B) is a solution because 

/■oo 

Axa{B) + xa{B)A = / e-'\AB + BA)e-'^dt 

Jo 

Conversely, suppose that AX + = B, then 

/■oo 

Xa{B)= e-'^{AX + XA)e-'^dt 
Jo 

{e-*^Xe-'^)dt = X , 



Jo dt^ 



'o dt 

and thus any solution X is equal to xa(-B). Hence the solution is unique. □ 
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